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D , Abstract 



We show that standard next-to-leading order (NLO) perturbative QCD analyses 
used to extract Ug from LEP data do not serve to disentangle the completely un- 
^ . known renormalization scheme (RS) invariant next-NLO (NNLO) and higher-order 

^ \ uncalculated corrections from those dependent on the renormalization scale in a 

predictable manner. The resulting quoted values of as{Mz) {^-Mg) with attendant 
scheme dependence uncertainties therefore reflect the ad hoc way in which they were 
extracted, rather than the actual values of these parameters. Choosing the scale so 
that the NLO coefficient vanishes (effective charge scheme) exposes the relative size 
of these unknown RS invariant higher order terms. They are seen to be sizeable for 
some of the LEP observables, indicating that they must be estimated if is to be 
determined reliably. This can be accomplished either from NNLO calculations, at 
present only available for the hadronic width of the {Rz), or non-perturbatively 
by writing the RS-invariant uncertainty in A in terms of the running dR{Q)/dlnQ 
of the observable R(Q) with energy using an effective charge formalism. The NNLO 
calculations for Rz and LEP data supplemented by lower energy PETRA data lead 
to Ai^ = 287 ± 100 MeV. We also discuss how the effective charge approach can 
be used to remove scale dependence from next-to-leading logarithm resummations 
of some LEP observables. 



1 Introduction 



The LEP collaborations now all have high-statistics data samples which enable them to 
make accurate measurements of a wide range of e^e~ QCD observables [1-9] - jet frac- 
tions, thrust distributions, the hadronic width of the energy-energy correlations, etc.. 
For most of these quantities non-perturbative effects such as hadronization corrections 
are expected to be reasonably small and next-to-leading order (NLO) calculations are 
available in renormalization group (RG) improved QCD perturbation theory [|lOi [TI| . For 



the hadronic width of the Z° alone we have a next-NLO (NNLO) calculation available 

By comparing such calculations with the data one hopes to be able to extract Aqcd, 
the fundamental SU(3) standard model parameter, or equivalently as, the MS scheme 
coupling constant. 

Unfortunately this intention is hampered by our rather limited current knowledge of 
what QCD actually predicts for any of these quantities due to the dependence of pertur- 
bative calculations at NLO on the unphysical renormalization scale fi, and more generally 
on the renormalization scheme (RS) in higher-orders. In a complete all-orders calculation 
of any quantity the /i-dependence would cancel and so it is an arbitary parameter, but 
our reliance on a truncated expansion in 05 in which such a cancellation no longer quite 
occurs renders predictions ambiguous. 

Since it was first clearly formulated in ref. this 'scheme dependence problem' 
has been the subject of an extensive literature which has yet to produce consensus. The 
response of theorists to the dependence of NLO predictions on /i has mainly been to 
advocate - in the face of the central fact that it is actually an arbitrary unphysical quantity 
- that some particular choice of the scale is the 'correct' or 'best' one, thereby removing 
the ambiguity. Any such choice must appeal to some criterion outside of the familiar, 
well-defined business of calculating Feynman diagrams and no such appeal has so far 
been able to command universal support. Confronted with such confusion amongst their 
colleagues, experimentalists have tended to adopt pragmatic compromises that 'mix and 
match' different proposed solutions. Specific reservations about the way this has been 
done at LEP will be expressed in section 2, but it is clear that any approach to the 
problem that increases consistency amongst observables by expanding error bars to include 
a 'theoretical uncertainty' (essentially obtained by averaging across camps of theorists) 
does so by sacrificing unambiguous clarity over what is being measured, and limits the 
extent to which QCD could be investigated in detail. 

In some sense, all previously proposed approaches have been an attempt to second- 
guess the uncalculated higher-order terms in the perturbation series, in that they say 
this choice of scale and a NLO calculation will give a prediction as if we knew the full 
series - and it is likely to be a reasonable, perhaps good, prediction at that. On the 
other hand, our position is essentially that the higher-order terms are simply unknown 
and that we therefore should restrict attention to the more modest aim of deciding what 
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it is actually possible to learn about the uncalculated terms from the information at our 
disposal. A widespread belief already exists that the severity of the scale dependence for 
any particular observable is somehow related to the size of its corrections in higher-orders, 
although this overlooks the problem that the size of the unknown terms is itself a scheme 
dependent notion - there is after all always a scheme in which there are no corrections 
whatsoever. Put another way, the total uncalculated corrections to a /i-dependent NLO 
prediction must also depend on the scale so that their sum does not. Superficially, we have 
therefore posed ourselves an unphysical (and hence uninteresting) question. However this 
/i-dependence would only be a fundamental problem if it were unknown, unpredictable 
and uncontrollable, when in fact (as we shall show) the higher-orders can be split into 
a predictable contribution incorporating all the awkward /i-dependence and a remaining 
piece containing all the genuinely unknown aspects. 

This is really just a consequence of the fact that the /i-dependence in a NLO calcu- 
lation has a universal form, and so a minimal version of our claim would be that we are 
suggesting a simple, convenient means of comparing theory and data in which the scale 
dependence is trivial. The formalism to be proposed has the advantages of being derived 
non-perturbatively and as far as possible entirely in terms of RS-invariant quantities. In 
particular, is avoided in favour of A;^, the dimensional transmutation parameter of 
QCD, as the free parameter in the theory since the former is just a theoretical construct 
without any direct experimental significance and is indeed strictly a different quantity in 
NNLO from what it was in NLO. It remains true that A is RS-dependent, but with the 
difference that the values in different schemes are exactly related by only a NLO calcula- 
tion. This A is to be as directly related as possible to R{Q), any generic LEP observable 
dependent on a single dimensionful scale Q (i.e. a/s, the e^e~ centre of mass energy). 
Because this R{Q) will play a role analogous to as, we shall refer to this approach as the 
'effective charge formalism'; it was originally introduced by Grunberg ||T^ and has also 
been discussed by Dhar and Gupta [|l^ . Unfortunately, there has been a widespread view 
that Grunberg's approach constitutes no more than a particular choice of RS. Indeed the 
insistence on using the RS-invariant R{Q) rather than the RS-dependent in the role 
of coupling means that the formalism has much in common with the fastest apparent 
convergence (FAG) criterion for choosing the scheme, by which the RS is chosen such 
that the higher corrections vanish. This perturbative property is a nice one, often simpli- 
fying calculations, but in this paper it is entirely secondary to the clarifying benefits of 
non-perturbatively identifying each observable as a coupling (from which it follows triv- 
ially). We cannot overemphasize that the effective charge (EG) formalism is both derived 
non-perturbatively and is more general than the adoption of a particular scheme, even if 
pieces of that formalism can be interpreted in terms of the FAG scheme. 

Because it should provide a way of measuring A without any scale ambiguity whatso- 
ever, thereby saving the formalism from merely being a useful framework for displaying 
data, as important an interpretation of the formalism derives from the way in which R{Q) 
runs as a renormalisation group improved coupling with Q. This evolution is described 
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by a function p{R) which is simultaneously both the basis for a convenient RS-invariant 
measure of the importance of the higher-orders and the /3-function corresponding to the 
FAC scheme. Thus the effective charge /3-function p{R) can, at least in principle, be 
experimentally measured, so that by detailed investigations of the running of R{Q) with 
Q one can determine even in the absence of NNLO perturbative calculations. The 
formalism applied in this way is non-perturbative and transcends the use of any particu- 
lar scheme, with the only barrier to an accurate determination of A;^ being the purely 
experimental uncertainty. Of course at present we do not have high precision information 
on the running of these LEP observables, and the running can only be investigated using 
data obtained at lower energies with different detectors and machines. In an ideal world 
LEP would take data at energies below y/s = Mz and study the running of the observ- 
ables in detail. As we shall argue, this would enable reasonably precise determinations 
of AjYg, and more importantly would also indicate whether NNLO corrections by them- 
selves would be sufficient, helping to establish whether the extremely arduous calculations 
involved are worth embarking on [|l^ . 



Even in the absence of information on NNLO corrections and the detailed running 
of the observables with energy, however, there is still information to be obtained with 
the effective charge formalism, given only NLO calculations. One can exhibit the relative 
size of the uncalculated higher-order corrections for different quantities. In the effective 
charge language these are related to how the energy dependence of the quantity at the 
experimental energy differs from its asymptotic dependence as Q —>■ oo. If this difference 
cannot be neglected, then NNLO calculations or measurements of the running are neces- 
sary before A can be determined. This can be decided in the effective charge formalism 
given only a NLO calculation and data at one energy. 

Our proposal is thus not a solution to the scheme dependence problem in the conven- 
tional sense. For a start, sufficiently much can be learnt about QCD independently of 
the absolute value of Ajjg for the problem itself to be something of a distraction; refining 
the experimental data on the relative size of higher-order corrections for different observ- 
ables and providing explanations for the aberrant ones is a serious challange in its own 
right. Furthermore, the energy-dependence data already allows us to make tentative ab- 
solute statements in a way complementary to perturbation theory. The scale dependence 
problem is not so much solved as completely outflanked. 

A preliminary discussion of the effective charge formalism applied in this way has 
appeared in references [17] and [18]. 

The plan of this paper is as follows. In section 2 we shall review the RS dependence 
problem and outline the traditional approach to as extraction. In section 3 we shall 
introduce in detail the effective charge formalism sketched above and stress its advantages 
over the conventional approach. We shall attempt to determine Aj^g using the rather 
limited information on NNLO corrections and energy dependence of observables currently 
available. We shall also discuss an attempt to resum higher-order terms in the effective 
charge beta-function using the knowledge of leading and next-to-leading logarithms in 
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the jet resolution variable yc for jet rates in the or Durham algorithm where 
the perturbation series has an exponentiated structure [1^, |2l|. Such a next-to- leading 



logarithm (NLL) resummation of the effective charge /5-function is without ambiguity, and 
in the effective charge language the Uc dependence can be extracted from the experimental 
measurements and compared with that of the resummed /3-function, indicating whether 
or not the NLL approximation is adequate. Recent attempts to compare resummed 
perturbation series with the LEP data require use of an ad hoc matching procedure 
|l|, 20| to include the exact NLO results, the additional problem of scale dependence still 
remaining. 

In section 4 we give our conclusions. 



2 Review of the Scheme Dependence Problem 



A. Parametrizing RS Dependence 



Consider a generic dimensionless LEP observable R{Q), where Q denotes the single 
dimensionful scale on which it depends, typically the e~^e~ centre of mass energy y/s. In 
RG improved perturbation theory we can write without loss of generality, 

R{Q) =a + ria'^ + r2a^ + .... (2.1) 

where a denotes the RG improved coupling a = ag/iT. Notice that by dividing any 
observable depending only on a single dimensionful scale by its possibly dimensionful 
tree-level perturbative coefficient (which will be RS invariant), and raising to a suitable 
power, we can always obtain a series of the form (|2.1| ) representing a dimensionless R{Q). 

The coupling a and series coefficients depend on the RS employed. Using the 
compact notation and conventions introduced by Stevenson ||13[ a is specified by a P- 
function equation 

^ = -a^{l + ca + + ...) = -/3(a), (2.2) 

where r = &lnp with /i the renormalisation scale and A the dimensional transmuta- 
tion mass parameter of QCD, here differing from the traditional definition by a factor 
(2c/6)-^/^ h and c are RS invariants dependent only on the number of quark flavours 
and the number of colours Nc\ for A^^c = 3 we have h = (33 — 2Nf)/6 and c = 
(153 — 19A^/)/2(33 — ^Nj). In all our LEP determinations we take Nf = 5 and as- 
sume massless quarks. For massive quarks the scheme dependence discussion will also go 
through in any mass-independent RS, i.e. one in which the coefficients C2, C3 . . . do not 
depend on the fermion masses . 
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The higher coefficients C2, C3, ... are RS-dependent. Indeed Stevenson [|T3|, extending 
the work of Stueckelberg and Petermann p3| , has shown that one may consistently use the 



parameters r, C2, C3, . . . to label the renormalization scheme. In the conventional approach 
when retaining terms up to and including r^a"^^ in equation one truncates the 

/^-function of equation ( |2.2| ) retaining terms up to and including c„a"'*'^. On integrating 
up the truncated equation (|2?^ ) one can define a^^^r, C2, . . . , c„) and correspondingly one 
finds for consistency ri(r),r2(r, C2),r3(r, C2, C3), . . .,r„(r, C2, C3, . . . , c„). In this way the n^^- 
order truncated approximant is also labelled by the scheme variables, R^^\t, C2, . . . , c„). 
Of course when summed to all orders this dependence must cancel and a formally RS- 
independent sum be obtained. 

The formal consistency of perturbation theory further ensures that 

i?(")(r', 4, . . . , 4) - i?(")(r, C2, . . . , c„) = ka^"-' + 0(a"+^), (2.3) 

so that differences between results in two different RS's are formally effects one order 
higher in perturbation theory. Of course k depends on {r, C2, . . . , c„} and {r', c'2, . . . , c'^} 
and may well not be a small coefficient. 



B. The RS dependence of R^^\t). 



Before discussing the general situation further let us consider the simplest n = 1, 
NLO case. This is the accuracy to which all but one of the LEP observables have been 
calculated. 

To NLO we have 

rW(t) = a«(^) + ri(r)(a(i)(r))2 (2.4) 
where a^^^(r) is obtained by integrating up the NLO truncation of equation (|2.2| ) 

^ = -a^n + ca). (2.5) 
dr 

With the boundary condition a^^\0) = 00 one obtains 

where we define the function F for later use, so that a^^^r) = F~^{t). 

How does ri(r) depend on r explicitly? To see this, consider two different RS's, RS 
and RS'. The connection between the corresponding renormalised couplings a and a' will 
then be such that 

/3(a) = ^P'{a) (2.7) 
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and we can define 

a' = a(l + z/ia + z/2a^ + . . .)• (2.8) 

Inserting ( p.8| ) into the series for R{a') and R{a) in tlie two RS's and equating coefficients 
one finds 

ri = z/i+r;. (2.9) 
By integrating up f3{a) and (3' {a') as in ( p.6|) and taking the difference we find 

t-t' = i + clnf-^) +0(a) 

_1 _c In f^^^ +0(a') (2.10) 

where the 0(a) and O(a') terms reflect contributions beyond NLO in the /3-functions. 
Using equation (|2.8|) and equating coefficients of corresponding powers of a gives 



t-t' = vi. (2.11) 
Ehminating Vi between equations (2]^) and ( 2.11 ) we flnally arrive at 

r[-ri = T' -T. (2.12) 

This implies that 

ri(r) =r + ri(0). (2.13) 



It is useful to identify the RS-invariant combination [|1^ 

Po = r-ri(r). (2.14) 
Since R{Q) is a function of a single dimensionful scale Q it follows that 

Po(Q)=r-ri(r)^Mn|, (2.15) 

with po{Q) and A RS-invariants. A is dependent only on the particular QCD observable 
R. This strongly suggests that these quantities should have a physical signiflcance as 
opposed to RS-dependent quantities such as ri(r) and a(r) which depend on unphysical 
parameters. As we shall show in section 3 this is indeed the case, po{Q) and A are 
connected with the asjTiiptotic {Q oo) dependence of R{Q) on Q with 

RiQf-"^^, = 77^, (2.16) 
poiQ) 6 In I 

and so Po{Q), or equivalently A could in principle be directly measured given unlimited 
experimental energies. 
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This observable-dependent 'physical' quantity A is directly connected to Ars the uni- 
versal dimensional transmutation parameter of QCD which we are attempting to deter- 
mine. To see this, recall that 

r = 61n^ (2.17) 

Ars 

is the variable which specifies the RS in NLO. Notice it is not sufficient to specify that 
one has chosen a given renormalisation scale n in order to specify the RS even at NLO. 
The renormalisation scheme is specified by r which involves Ajis as well. 

Since Pq{Q) is an invariant we see that for two different RS's, RS and -RS" 

Mn J^-rf^(/i) =6 In J^-rf^'(/i), (2.18) 
Ars Ars' 

where rf'^(/i) denotes the NLO correction evaluated in the renormalisation scheme RS 
and rf"^(yu) = ri(r) with r = blnfi/XRs- From equation ( |2.15| ) with fj, = Q we find 
directly that 

A = AH5exp(rf (/i = Q)/6) (2.19) 

and so given a NLO calculation in some RS the invariant A is exactly related to the 
dimensional transmutation parameter Ars- Notice that rf'^(yU = Q) is a Q-independent 
quantity. 

Rearranging equation ( 2.18|) yields 



Ars> = A^sexp[(rf^(/i) - {^^))/b] (2.20) 

The exponent is a universal /i-independent constant which exactly relates the A's in the 
two schemes, given NLO calculations in these two RS's. The implication of this important 



result [13, 2^| is that it does not matter which A^5 we choose to try to extract from the 
data. If the RS is mass-independent then Ars will be universal and can be exactly 
related to Ars' in any other scheme by a universal factor given exactly if we have a NLO 
calculation for any observable in both RS's. Thus there is no non-trivial residual scheme 
dependence implied in our convenient choice of A;^ as the fundamental parameter we 
wish to extract. 

As an illustration consider the minimal subtraction [MS) renormalisation procedure 
and the MS or modified minimal subtraction procedure where the ln(4'7r) — 7^; terms 
present in dimensional regularization are subtracted off (7^; = 0.5772 ... is Euler's con- 



stant). For any observable, and independent of 



(/i) - rr (/i) = ^(ln(4vr) - je). (2.21) 



So using equation ( |2.19| ) we have 



Ams = AMsexp(i^^^^^^^^) = 2.66Ams. (2.22) 



7 



This relation is independent of Nf or Nc. For the MS and MS procedures the higher 
order RS-dependent /9-function coefficients are identical, c^^^ = c^^^ {k > 2). The only 
difference is in the subtraction procedure and hence in the definition of the renormalisation 
scale fi. A choice of scale using the MS procedure corresponds, therefore, to the same 
renormalization scheme as use of the MS procedure with scale 2.66yU. 

We are finally in a position to exhibit the explicit dependence of R^^\t) on r. In fact 
it is easier to equivalently consider the dependence on a^^^r). Using equation (|2.15|) to 



write ri(r) in terms of po ^^nd r, and using for a^^^r) in terms of r we find 

J?(i)(r) = a«(r) + (a«(r))2F(a«(r)) - po(Q)(a«M)'- (2.23) 

When plotted versus a(r), R^^^ has the generic approximately inverted-parabolic shape 
shown in Figure 1, provided that Pq > 0. From ( |2.15|) this condition will automatically 
hold provided that Q > A. 



C. RS dependence of r„(r, C2, . . . , c„) 



Before we discuss this NLO scheme dependence in more detail let us complete the 
picture by discussing the explicit scheme dependence of the higher coefficients r-2(r, C2), 
i^3i'T,C2,C3), ... as well. Integrating up the /5-function truncated at n*'^ order with an 
infinite constant of integration (related to the definition of A, see section 3) we obtain, 
analogous to equation (Of) , 



1 



+ cln 



ca 



(n) 



+ 



dx 







X] 



(1 + CX + C2X^ + CsX^ + . . . + CnX^] 



1 



X 



+ 



1 



[l + cx) 



[2.2A) 



This transcendental equation can then be solved explicitly for a(")(r, C2 



To obtain the explicit RS dependence of the r, it is convenient to identify a special 

J, in which ri 



RS, the effective charge (EC) scheme [14 
Then R^'^^ 



r2 = rs = . . . = r„ = 0. 

a^"^ in this scheme, and the coupling constant is the physical observable. 
This choice of RS will correspond to particular values of the scheme parameters, tec = 

The fact that tec = Po follows directly from equation (|2.15|) , 'ri(r) = r — po then 
ri{TEc) = Tec — Po = =^ tec = Po- From equation ( p^) we have that for two RS's, 
barred and unbarred, 

da 



If the barred RS is chosen to be the EC scheme then 



P{a) = p(a) 



a^(l + ca + p2a^ + . . . + p^a^ + . . .) 



:2.25) 



(2.26) 
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with a = R. Then ( p.25| ) gives 



p{R) = (3{a{R)) — , 



(2.27) 



where a{R) is the inverted perturbation series. By expanding out the right-hand side as 
a power series in R and equating coefficients we obtain 



(2.28) 



P2 = C2 + r2 - ric - 

Ps = C3 + 2r3 - 4rir2 - 2rip2 - rjc + 2rl 



These objects are Q-independent RS invariants constructed from the and Cj in any 
RS. For instance given NLO and NNLO calculations of rf^^{Q),r2^'^{Q) and since Cg^*^ is 



known 27, 28 , one can construct 



P2 = cf^ + rf^_rf^5c-(rf^)2. 



(2.29) 



As we shall see in the next section just like the NLO invariant po, the NNLO and higher 
P2,P3, ■ ■ ■ invariants have a physical significance, whereas quantities such as rf^^{Q),c^^ 
are intermediate RS-dependent quantities which should be eventually eliminated in favour 
of RS invariants. 

Having obtained p2, P3 from NNLO, and higher order calculations in any convenient 
RS, we can exhibit the explicit r, C2, C3, . . . dependence of the perturbative coefficients by 
rearranging equation ( p.28|) . We have that 

n{r) = (r-po), 
r2{T,C2) = (r - po)^ + c(r - po) + (P2 - C2), 

r3(r,C2,C3) = (r-po)' + |c(r-po)' + (3p2-2c2)(r-po) + ^(P3-C3) (2.30) 



The result for r„(r, C2, . . . , c„) is a polynomial in (r — po) of degree n with coefficients 
involving p„, p„_i, . . . , c and C2, C3, . . . , c„, such that r„(po, P2, • • • , Pn) = 0. 



D. NLO extraction of ^^(M^). 



Let us now return to Figure 1 and discuss how one might attempt to use it to extract 
■^JIs cns{Mz) in the MS scheme. 

Notice ffist that the curve R^^^ (r) of equation ( 2.2'S ) is a universal function of r and 
Po (for given fixed Nj), where the value of the invariant po depends on the particular 
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observable R. po is of course directly related to given a NLO calculation of rf^^{Q), 
with 

Po = b\nJ--rY^{Q). (2.31) 

i?(i)(r) has a maximum at a(r) ~ 1/po (t — Po) where R^^\po) ~ 1/po- The full-width 
at half-maximum in a(r) is approximately \/2/po. -R*-^^ vanishes at a = (r = cxd) and at 
a(r) ~ 2/po (r ~ po/2). For comparison Figure 1 is drawn with po = 25 and A^j = 5. 

Let us suppose that the horizontal dashed line in Figure 1 represents the measured 
experimental data. If po is sufficiently small that the maximum R^^^ ~ 1/po lies above the 
data then the curve will cut the data at two scales ti,T2 (as illustrated). Conversely, if 
Po is made larger the curve will be below the data line. Thus an infinite set of r, po pairs 
fit the data perfectly. If we wish to measure po (and hence A^g^) we must specify r. This 
is the NLO scheme dependence problem. At NNLO one would have a surface R^'^\t^C2) 
and to extract po one would need to specify r and C2, and so on with one extra unphysical 
parameter for each order in perturbation theory. At least in NLO for a given value of po 
(A;^) there is a maximum possible R^^\ 

In NNLO and higher one can show that for a given value of po (A]g^) there exists 
a choice of r, C2, . . ., c„ such that i?*^"^ has any desired positive value. Thus, for any ^jj^ 
we can choose a sequence of schemes such that R^'^'^ = R^^^ = . . . = R^^'> = Rexp, the 
experimentally measured value. 

Various "solutions" of the scheme dependence problem, i.e. motivations for particular 
choices of r, C2, . . . , Cn have been proposed. 



(1) Principle of Minimal Sensitivity (PMS) ||T3[| 



The idea here is that since the exact (all-orders) result is independent of r,C2,C3,. . . c„ 
one should choose the n'^^-order approximation i?*^"')(r, C2, . . . , Cc) to mimic this property 
and be as insensitive as possible to the chosen value of the unphysical parameters. That 
is one arranges that 



dr 



QR{n) 



C2=C2 







and the PMS scheme is specified by r, C2, . . . , c^. At NLO one has to solve 



dr 



0, 



and this corresponds to solving the transcendental equation 



-|- cln 



ca 
1 + ca. 



1 c 

2 1 + ca 



Po- 



(2.32) 



(2.33) 



(2.34) 
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The solution a(r) = a then gives 



(1) '^(1 + |ca) 
RpMS — /I I =\~- (2.35) 
(1 + ca) 

The stationary PMS point of R^^\t) is the maximum at r = r ~ po ^ c/2 (this is a 
somewhat better approximation than the cruder r ^ po given previously for the position 
of the maximum), po is adjusted so that RpMs — -^e 



(2) Effective Charge (EC) Scheme [|T|, p 



This corresponds to a choice of scheme such that ri = r2 = . . . = r„ = 0, hence 
is an effective charge. As we have seen the scheme parameters are then 
{po,p2,P3, • • • ,Pn}, where the pi are the RS-invariants in equation ( |2.28| ) . Once these 



have been computed up to pn from higher-order calculations in any technically-convenient 
RS (e.g. MS) they can be inserted in the integrated /9-function equation (2.24) which can 
then be solved with a*^") = R for given po. po can then be obtained by requiring R = Rexp- 

This particular approach is sometimes also referred to as the fastest apparent conver- 
gence criterion. In this language it appears rather artificial, but as we shall emphasise 
in section 3 its advantage is that the all-orders coupling constant and /^-function are 
experimentally observable physical quantities, allowing a non-perturbative approach to 
supplement measurement, and a physical definition of the uncertainty. 

It should be emphasised that EC and PMS predictions in NLO are very similar since, 
as we have seen, tec = Po, tpms — Po — c/2. It is also true that in NNLO EC and 
PMS remain close to each other [RD, pT|. For the PMS approach, however, the coupling 



constant and /3-function are unphysical quantities, and it is not clear even that their all- 
orders versions are defined, given the complex nature of the coupled equations which must 
be solved. 



(3) Physical Scale |P2 



According to this viewpoint the renormalization scale should be chosen to be close 
to the physical scale in the problem, p = Q. If predictions in the vicinity oi p = Q are 
strongly p-dependent then this is supposedly an indication that the perturbation series is 
intrinsically badly behaved. This viewpoint is usually justified by noting that perturbative 
coefficients in higher orders will be polynomials in Inp/Q, 

rn = Y.Kni{b\n^y. (2.36) 

1=0 ^ 
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To avoid unnecessarily large logarithms one should therefore set /i ~ Q. Whilst this is 
true as far as it goes, it tacitly assumes that the NLO RS dependence is completely given 
by the dependence on the renormalisation scale. In fact as we have seen r„(r, C2, . . . , c„) 
where r = hln^/KRs so the coefficients Kni in (|2.36| ) will depend on Krs- The meaning 
of n of course depends for instance on whether it is fj,{MS) or fi{MS). In contrast we 
have from (|2.30|) that 



l^, C2, 



1=0 



(2.37) 



So to avoid unnecessarily large terms we should clearly choose r ~ po? the effective charge 
scheme! To make contact with ( 2.36 ) notice that 



- Po) 

So we can write (|2.37| ) as 



Mn^ + rf (Q) 



6 In 



, 1^^, . . . , 



1=0 



bin 



fie 



Q 



(2.38) 



Q 



^2.39) 



In (|2.37| ) and (|2.39|) the coefficients Kni do not depend on the NLO RS choice, only on 
C2, C3, . . . and the RS invariants p2, P3, ■ ■ ., whereas in (|2.36| ) the Kni will depend on rf'^(Q) 
and so have a hidden dependence on the NLO RS choice which is customarily ignored in 
the usual 'physical scale' argument. Applying the argument to ( p.39| ) instead one would 



infer that one should set p ~ Qe 



'{Q)/b 



QAjis/A, which of course corresponds to the 



effective charge scheme, r = po- For particular cases the extra factor may be large. 

Recall that A is an RS-invariant characterising how the observable runs with Q asymp- 
totically. Different observables will have different A's. If it happens that Ars — A 
(rf'^(Q) ~ 0) then predictions will be stable and large logarithms avoided for fi c::^ Q, but 
no special physical significance should be accorded to this fact. It is always possible to 
modify the subtraction procedure so that Afjs = A, since A's in different RS's are related 
according to ( |2.2q ) . 

We conclude that a modified 'avoidance of unnecessarily large logarithms' physical 
scale argument which correctly labels the NLO RS dependence actually leads to the 
effective charge scheme, r = po- 



(4) Fitting p and A to the data 



For most of the LEP observables we have a differential distribution in a kinematical 
variable rather than a single experimental data value as in Figure 1. For instance a 
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thrust distribution in the thrust variable T, or a multijet rate in the jet resolution cut 
He. Denoting such a generic kinematical variable by A, we should really add an extra 
perpendicular axis to Figure 1 and consider the A dependence as well, R^^\t,\). The 
NLO coefficient will also depend on A and so for a given choice of ^jjs have po{\) 

from ( p^.31D . If is sufficiently large then the experimental data Rexp{^) will intersect 
R^^^T, A) for at least one r value for any A and so the data can be ffited perfectly with a 
suitable A-dependent RS parameter choice t(A). In general there will be two intersection 
points as in Figure 1 and so two functions ri(A),r2(A) will ffi the data perfectly. Since 
this can be done for any suitably large there is no unique best ffi. 

The OPAL analysis [|I], 0] uses a mixture of approaches (3) and (4). The data for each 
observable over a range of the corresponding kinematical variable is compared with the 
NLO prediction with the MS scale chosen as /x = Mz, and a best ffi for Ajjg performed. 
The resulting as{Mz) values are shown in Figure 2(a), reproduced from reference P]. 
There is evidently a considerable scatter in the as values obtained. 

One then performs a simultaneous best fit for /i and Aj^ for each observable over a 
range of the kinematical variables. This provides a second as{Mz) value. The quoted 
central value of as{Mz) is then taken to lie mid-way between these extremes (perhaps 
modified by hadronization) and the difference between them is taken as indicative of the 
size of the uncalculated higher-order corrections. By enlarging the 'scheme dependence 
error' in this way one obtains greater apparent consistency between different observables 
as shown in Figure 2(b). 

In our view one learns nothing useful from an analysis of this kind. Consider Figure 
2(a) and the considerable scatter of as values obtained from different observables ffited 
with II = Mz- Choosing /i = Mz means that from ( p.39|) the higher-order coefficients 
are given by 

n 

r„(/i = Mz) = Y.KUrf''{Mz)y (2.40) 

where the Kni depend on the invariants P2, Ps, ■ ■ ■ , Pn and on the NNLO and higher RS 
parameters C2, C3, . . . , c„. We have tabulated in Table 1 the rf^'^(M^) values for some of 
the observables at particular values of the kinematical parameters using the NLO matrix 
elements of reference 10. (For more details see section 3B). One sees that they are large 
and rather scattered, e.g for the jet rates rf^^{Mz) ~ 10 which results in rather small EC 
(or PMS) scales, pec ~ 0.07M^. As emphasized earlier no physical significance attaches 
to these small scales, since the scale itself is just an artifact of the unphysical subtraction 
procedure employed. Thus from ( p. 401 ) one may expect large and rather scattered higher 



order corrections for the different observables with p = Mz, which translates itself into 
the expectation of scattered values of as{Mz), which indeed are observed. 

If we instead choose p = Mze~'^^^^^^^^^ = p^c (the effective charge scale) for each 
observable then the higher order coefficients are given by 

Tnip = Pec) = (2.41) 
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where Kno does not depend on the NLO RS choice (r), only on the NNLO and higher RS 

parameters and the RS-invariants p2,p3, . . .. Any scatter now observed in the extracted 
^Ws '^^ (^siMz) for different observables will be attributable to these NNLO and higher 
RS invariants being large, and not to the scatter due to the NLO scheme dependence 
logarithms involving r^^{Mz), which are avoided by the choice n = Hec- As we shall 
discuss in Section 3 there is still scatter when this is done but it is somewhat reduced. 

With the choice /x = Mz predictable scatter due to avoidable already known NLO 
scheme dependence logarithms is superimposed on top of the interesting scatter due to 
the size of the NNLO and higher invariants which is telling us about NNLO and higher 
order uncalculated corrections. The choice n = Hec removes the predictable scatter and 
provides genuine information on the interesting NNLO and higher-order uncalculated 
corrections. 

The MS scales obtained by fitting to the data for an observable over some range of 
the associated kinematical parameters tend to be considerably smaller than Mz, as for 
the EC or PMS scale. Typically for values of the kinematical variable A well away from 
the 2-jet region where large logarithms will be important, po{X) and ri(A) do not vary 
strongly with A. This means that one cannot obtain best fits over such a range of A. It is 
only by including the 2-jet region that one can obtain stable best fits, but then fitting for 
a A-independent scale over a range where po(A) and ri(A) (and presumably uncalculated 
higher order corrections) are strongly A-dependent represents a rather severe compromise, 
and it is not obvious what one learns. 

In conclusion, we believe that the standard NLO analyses used to extract as from 
LEP data do not serve to disentangle the genuinely new and interesting NNLO and 
higher order uncalculated corrections from predictable higher order corrections connected 
with the choice of RS at NLO. The resulting quoted values of as{Mz) (or A.jjg) with 
attendant scheme dependence uncertainties therefore reflect the ad hoc way in which 
they were extracted, rather than the actual values of these parameters. In contrast, by 
choosing the effective charge scale for each observable one exposes the relative importance 
of uncalculated NNLO and higher corrections. An effective charge formalism which can 
be supplemented with non-perturbative information is described in the next section. 



3 The Effective Charge Formahsm 

A. The Q dependence of R{Q) 
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For the dimensionless QCD observable R{Q) we can define 

dR 



dlnQ 



^m- (3.1) 



dR/dlnQ and hence ^{R{Q)) are, at least in principle, experimentally observable quan- 
tities, although collider experiments are usually designed to make high-statistics observa- 
tions at a fixed energy Q rather than examining the detailed running of observables with 
energy. To make contact with the standard perturbative approach we note that ( p.l|) is 
the /3-function equation in the effective charge scheme where the coupling a = R{Q) and 



the /3-function is p{R) as defined in equation ( p.26| ) . So we have 

dR 
d\nQ 



e(i?(Q)) = -bpiR) ^ -hR\l+cR+p2R^+. . .)+e-^l'^R\Ko+K^R+. . .). (3.2) 



The effective charge /3-function p{R) may be regarded as a physical observable. As mea- 
sured from data it will include a resummed version of the (perhaps asymptotic) formal 
perturbative series exhibited in ( p.2| ) together with the non-perturbative terms involv- 
ing e~^/^ which are invisible in perturbation theory. (Here S,5,-ft'o?-^iv • - are observable- 
dependent constants which we shall not specify further [^.) Hence by comparing with 
data one can test how well the first few terms of the perturbative series serve to rep- 
resent the observed running of R{Q). That is, from measurements of RexpiQ) and 
dRexp{Q)/dlnQ, one can test the extent to which 

^ - -bRlpil + cRexp) (3.3) 

describes the data. A marked discrepancy would either indicate the importance of the 
NNLO invariant p2 and higher terms, or the relative importance of the non-perturbative 
contributions, or both. 

Given a collider program dedicated to the investigation of the Q-dependence of observ- 
ables the above tests would indeed be powerful and useful. The focus of present studies, 
however, is the accurate measurement of R at fixed Q, i.e. Q = Mz, and the comparison 
of these results with NLO matrix element calculations in order to extract Ajjg. To obtain 
R{Q) we clearly need to integrate equation (|3.1|) . The boundary condition will be the 
assumption of asymptotic freedom, that is R{Q) — as Q oo, which corresponds 
to the requirement that ^{R{Q)) < for Q > Qq, with Qq some suitably low energy. 
(Equivalently, p{R{Q)) > for Q > Qq, assuming 6 > or A^j < 33/2 for A^c = 3 QCD.) 
Any zero of ^(i?(Q)), ^(-R*) = 0, say, would correspond to R{Q) —^R*asQ^C!C and 
ultraviolet fixed point behaviour. 

Integrating equation ( |3.1| ) we obtain 

In = = / dx + (infinite constant) (3.4) 

A Jo t,{x) 
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where A is a finite constant which depends on the way the infinite constant is chosen. 
The infinite constant can be chosen to be 



{Infinite constant) 



dx 

■r]{x) 



(3.5) 



where ri{x) is any function which has the same x —>■ behaviour as ^{x). We know from 
equation (|3.2|) that ^{x) has the universal a; ^ behaviour C,{x) = —bp{x) ~ —bx'^{l+cx) 
so we choose ri{x) = —bx^{l + cx). Inserting this choice for r]{x) and rearranging equation 
Tl we find 



m4 

A 



dx 



R{Q) X^{1 + CX^ 



+ 



dx 



p{x] 



+ 



1 



x^fl + cx) 



(3.6) 



The first integral on the right-hand side is just F{R) where F is the function defined in 
equation (2.6), 



F{x) 

We define for later convenience 
Apo(Q) -= 



+ cln 



X 



cx 



1 + cx 



R{Q) 



dx 



1 



1 



(3.7) 



(3.^ 



p{x) x'^{l + cx) 

Notice that despite appearances the integrand of Apo is regular at x = 0. Then we have 



F(i?(g)) = Mn|-Apo(g). 
As Q oo, R{Q) — so that Apo and thus asymptotically 

Q 



F{R{Q)) ^b\n 



A 



with the constant of integration given by 

A= lim Qexp{-F{R{Q))/b). 



(3.9) 



(3.10) 



(3.11) 



Notice that in arriving at (3.9) we did not need to refer to perturbation theory, except 
to assume the asymptotic a; — > (Q ^ oo) behaviour ^{R{Q)) f» — 6i?^(l + cR). Despite 
this, the constant A obtained with the particular choice 7]{x) = —6x^(1 + cx) is precisely 
the A introduced in equation ( p.l5|) . That is 



A = A^exp(rf^(p = g)/6). 



(3.12) 
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We can see this immediately since if a denotes the couphng in the RS with jj, = Q using 
the MS subtraction procedure in NLO, and with higher order /3-function coefficients zero 
p6| we have 

h\i^S— = F{a) (3.13) 
giving a well-defined all-orders coupling. Further 

R = a + rY^{Q)a^ + .... (3.14) 
a — > as Q — > oo, so asymptotically as Q — > C)0 we have 

F{R)^ F{a)-rf^{Q) + ... (3.15) 



where the ellipsis denotes terms which vanish as Q ^ oo. Inserting ( p.l5|) into ( |3.11|) 
and using ( 3.131) we find (|3.12|) as Q ^ oo. 



So finally we have 



F{R{Q)) = Mn|-Apo(g) 

Q .JTSr 



&ln^-rf^(Q)-Apo(Q) (3.16) 
PoiQ)~ApoiQ). 



We could, of course, have written down the result of equation ( p.l6|) perturbatively in 



the EC scheme at once by simply using the integrated beta-function equation of equation 
(2.24) with r = po, = p(") = l + cx + p2x'^ + . . . + pnX"", and a^") = However, our 
purpose here has been to stress that equation (|3.16| ) holds beyond perturbation theory 
for the measured observable R{Q) and Apo(Q) constructed from the measured running 
of R{Q), dR{Q)/d\nQ = —hp{R{Q)). That is, we can write a non-perturbative closed 
expression exactly relating the universal constant Ajjg to observable quantites. 

Since Apo(Qo) involves an integral in R{Q) between R{Qo) and 0, to actually measure 
it would require knowlege of R{Q) (equivalently dR/dlnQ) on the full range [Qo,oo). 
We of course only know R{Q) on some finite energy range, so if we want to obtain 
from equation ( p. 161 ) and measurements of R{Q) there will always be an uncertainty 



related to our lack of knowledge about the behaviour of R{Q) beyond the highest energy 
reached, or correspondingly related to our lack of knowledge about p{R) for R < R{Qo), 
with Qo the highest energy reached. As we shall discuss in sub-section 3C we can use 
measurements of the running of R{Q) to determine p{R) in the vicinity of i? = R{Qo), 
and perturbative QCD calculations to determine p(-R) in the vicinity of i? = 0. Putting 
these two pieces of information together we can make unambiguous statements about the 
validity of perturbation theory and the nature of the function p{R), which can assist our 
attempts to determine Ajjg. 
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Notice the fundamental significance of the NLO perturbative coefficient rf^^{Q) in all 
of this. Knowledge of the full behaviour of R{Q) at large Q would allow one to extract 
A from equation , but this by itself is observable- dependent and so one would 

not test QCD. A test is only possible if one also knows rf^^{Q) (from Feynman diagram 
calculations) and can then obtain the universal Aj^ from equation ( p.l2 ) . 



The universal Ajrq of course depends on the number of active quark flavours, A. 



MS ' 



the RS invariants b and c occurring in equation ( |3.16| ) also involve Nf. Since Aj^g in 



equation ( p.l6 ) involves the Q ^ oo behaviour of R{Q), one might wonder what value 
for Nf should be taken. With three generations of quarks, Nf = 6 presumably represents 
the asymptotic number of active flavours and so, introducing the notation b{Nf), c{Nf), 
the function ri{x) in equation ( p. 51 ) should be chosen as ri{x) = — 6(6)x^(l + c(6)x). 
So flnally ( p.l6|) with the explicit Nf dependence exhibited becomes 



{r{Q) 
6(6) In ^) 



+ c(6) In 



c{6)R{Q) 
l + c{6)R{Q) 



(3.17) 



MS 



A^ 

MS 



(Q, Nf = 6) 



rR{Q) 



dx 



m 

^{x) ' x2(l + c(6)x) 



+ 



Obviously the strict asymptotic Q ^ oo behaviour of R{Q) and dR/dlnQ = C,{R{Q)) is 
an idealized concept since we could never actually measure it. For instance, it presumably 
makes no sense to consider QCD in isolation beyond the GUT energy scale. Rather than 
using ( |3.17| ) for LEP observables where we have Nf = 5 active flavours, it makes more 
sense to consider the equation for R{Q) in a world with A'^^ = 5 active flavours, replacing 
6(6), c(6), and A^ in ( |3.17| ) by 6(5), c(5), and A^. The decoupling theorem 
means that we can consider an Nf = 5 version of QCD for energies below top threshold. 
As we shall discuss we shall be interested in how the observed Q-evolution of R{Q) at 
LEP energies compares with its asymptotic Q oo evolution in a hypothetical world with 
Nf = 5 flavours. This will decide how accurately we can determine the universal constant 
A^i^. To have any chance of determining A^^ we would require measurements around 
the top threshold and at higher energies, which we do not currently possess. 

A further subtlety connected with the effective charge formalism concerns the fact 
that one can always deflne effective charges other than a = R. More generally one could 
deflne /(a) = R where /(a) = a + /2a^ + /ao^ + . . . is any analytic function of a. Chyla has 
suggested that this arbitrariness is just the scheme dependence problem in disguise . 
In our view this is incorrect (see also Grunberg's response [§5). In fact the arbitrariness 
in deflning the effective charge is completely equivalent to redeflning the observable to 
be i? = f~^{R), maintaining the standard deflnition of effective charge a = R. One can 
of course always consider R as the measured observable rather than R, and then NLO 
truncation of the perturbation series for R using some RS will give different results from 
those obtained with R. Such a redeflnition of the observable would only be useful if some 
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information on the uncalculated higher order corrections, i.e. knowledge of the function 
p{R), were available to inform the choice. 
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B. EC formalism in NLO - the Apo plot. 



Let us suppose that we have measurements for a number of LEP observables at a fixed 
energy {Q = Mz), and knowledge of their NLO QCD perturbative corrections r^^{Q). 
What can we learn? 

Recall from equation ( p.l6|) that 



F{R{Q)) = Mn^-rf^(Q)-Apo(Q) 

Ws 

= po{Q)-Apo{Q), (3.18) 
where b, c and rf^^{Q) are evaluated with Nf = 5. As Q ^ oo, Apo(Q) and 

F{Rm^PoiQ). (3.19) 
At sufficiently large Q, Apo ^ po and ( |3.19|) may be used to obtain Aj^. We have that 

Ag^ = Qexp[-(F(i?(Q)) +rp(Q) + Apo(Q))/6]. (3.20) 

Given only a NLO perturbative calculation we have no information on Apo((5). From 
the measured observable Rexp and the NLO perturbative coefficient we can then extract 

exp(Apo(g)/&)Ag^ = gexp[-(F(i?,.,(Q))+rp(g))/6]. (3.21) 

To relate this to the discussions of section 2, the quantity on the left of equation (|3.21|) 
is just the value of A;^ which would be extracted from the data at NLO choosing p = 
Mz exp{—rf^ /b) = Pec, the effective charge scale. If Apo(g) <C 1, then this NLO 
estimate will be accurate. Indeed the fractional error in A^g^ will be 



6K 




. /Apo(g)\ 


A 




1 exp(^ 6 ) 



(3.22) 



The size of Apo(M^) will depend on the QCD observable. If the A^^ values obtained 



MS 

from equation ( p.21D exhibit significant scatter, then this indicates that Apq{Q) is 
not negligible for at least some of the observables. Recall that Apo(g) is defined non- 
perturbatively in terms of the running of R{Q) by equation (3.8) where dR/dlnQ = 
^(-R(g)) = —bp{R{Q)). If g is sufficiently large that dR/d\nQ has its asymptotic {Q 
oo) g-dependence, that is if p(-R) ~ R^^l+cR), then Apo(g) ~ 0. Thus signiffcant scatter 

~ ('5'} 

in the A^^ extracted at NLO with the effective charge scale unambiguously indicates that 
the current experimental Q is not yet large enough for all observables to be evolving with 
g asymptotically. If we want to reliably determine A^^ additional information on the 



'sub-asymptotic' effects is required. 
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The presence of these effects is a physical fact which can be directly observed from 
the Q-dependence of data, and cannot be remedied by changing the unphysical rcnormal- 
ization scheme. We stress once again how different this is from the scatter obtained with 
H = Mz which will be partly due to different values of r^^{Mz) and hence which can be 
modified by changing the unphysical subtraction procedure. 

To make this distinction more precise, let us define A(NLO, ri) to be the value 

extracted at NLO using an RS with rf ^(At) = n. The NLO coefficient ri may be taken 
to label the RS at NLO. Using a to denote the couphng constant in this RS we have 

R^a + na'^, (3.23) 

and 

bin ^ ^ = F(a). (3.24) 

AiNLO,ri) ^ ^ ^ ^ 

Then using equations (3. 20), (3. 23), and (3.24) it is straightforward to show that 

A{NLO,n) = exp[(AF(ri,i?) + Apo)/&]Aj^ 

= exp(AF(ri,R)/b)A(NLO,0), (3.25) 

where 

AF(ri, R) ^ FiR) - F (zl±^^±lF^ + r„ (3.26) 
AF(0, R) = 0. 

A{NLO, 0) is just the extracted choosing n = hec {'Ti = 0), given by equation (3.21). 

Prom equation (3.25) we see that the NLO scatter in is factorized into two 

components. The first is an already known ri-dependent contribution involving AF(ri, R) 
which may be obtained exactly from (3.26) with R = Rexp, and does not involve the 
unknown NNLO and higher RS invariants. The second component, Apo ~ P2R, does 
not depend on ri but does involve the unknown higher invariants, and represents the 
irreducible uncertainty in determining Aj^g at NLO. The uninformative ri-dependent 
part of the scatter can therefore be completely removed by choosing /i — /iec, which sets 
AF = 0. The scatter is then completely given by Apo- 

Rather than focusing on the scatter in A;^, we prefer for presentational reasons to 
concentrate on the relative size of Apo- We can define for each observable 



Ap^^iQ) = bln^- rf^ (Q) - F{ReM)- (3-27) 

MS 



Given a NLO perturbative calculation oirf^^{Q) and experimental measurements oi Rexp 
we can then measure Apo from the data up to b\nA^^ which should be a universal 
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constant. Changing merely translates the Apo obtained up or down by the same 
amount for each observable, so we can choose an arbitrary reference value for A^g^. For two 
different observables A and -B, the difference Apo'^^^A — ^Pff^^B can then be measured 
absolutely and so the relative size of Apo may be investigated. If Apo ~ for each 
observable, corresponding to small sub-asymptotic effects, then the Apo^^ points should 
lie on a horizontal straight line with some small scatter. Substantial deviations from a 
horizontal straight line therefore unambiguously indicate the presence of sizeable sub- 
asymptotic effects. The scatter in Apg^^ is of course entirely equivalent to the scatter in 
the A^r/j- extracted at NLO with u = llec^ these correspond to the reference A^r^ values 
for which Apo^^ = 0. We shall return to these A-^ values later. 

If one has a NNLO perturbative calculation or measurement of Rexp at more than one 
energy then, as we shall discuss in the next subsection, one can attempt to estimate Apo 
directly. One can then see if for all observables a single value of A^^^ brings Apn^^ into 
agreement with the directly estimated Apo. 

We now turn to extracting Apo^^ from LEP data. The observables we shall use are 
the ones appearing in table 1. R2 and -R3 denote the jet fractions defined by R2 = 
o'2{,yc) / CTtot and R3 = CT3{yc) / o'tot, where <72{yc) and (Xad/c) are the cross-sections for e~^e~ — > 
2, 3 resolved jets respectively, with i/c the jet resolution cut and atot the total hadronic 
cross-section. The EO and D labels denote different recombination algorithms used to 
cluster hadrons together. Details of these algorithms and results for the NLO perturbative 
coefficients rf^^{Q) can be found in reference []TT| for the EO algorithm, and for the D 



(Durham or kr) algorithm in reference |T9 



T denotes the thrust variable, and x the energy-energy correlation (EEC) angle. We 
used recent OPAL data on the differential distributions in these quantities |38|. The 
asymmetry in the energy-energy correlation (AEEC) obtained by subtracting the EEC 
measured at x and (180° — x) '^as also considered. Details of the definitions of these 
quantities and their NLO perturbative coefficients are contained in reference |0 . 

From the measured hadronic width of the (Thad) and the leptonic decay width 
(Tiep) one can obtain Rz = Thadf^iep = _(19.97 ± 0.03) (1 + 6qcd), where 6qcd has a 
perturbation series of the form (2.1) with rf'^^{Mz) = 1.41. The electroweak contribution 
(19.97 ±0.03) is from reference [^. This coefficient is for massless quarks, if one includes 
heavy quark masses [|^, one finds Sqcd = (1.05a + 0.9a^) instead in the MS scheme 
with fi = Mz- Within the errors the values of Apg^^ obtained are are consistent, so 
we shall use massless QCD. For the experimental value we take the 1993 LEP average 
Rz = 20.763 ±0.049 §. 

As noted earlier the reference value of A^^ assumed for Apn^^ is irrelevant since we are 

Ivl o 

concentrating here on the scatter of Apg^^ values. We shall choose A^^ = 110 MeV which 
corresponds to the central value obtained from a non-perturbative lattice analysis of the 



IP-IS splitting in the charmonium system [^. Whilst at present limited by the use of a 



quenched approximation, such calculations could in future provide reliable estimates of A. 
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In Figure 3 we show the Ap^^^ values obtained for the LEP observables discussed above. 
We have taken for each observable the particular values of the kinematical variables given 
in Table 1. Of course for each observable we actually have a Apg^^ distribution in the 
associated kinematical variable. In Figures 4-7 we give the Apg^^ distributions for the 
jet rates R2{Eo), R^^Eq), i?2(-D), Rsi^D) as functions of the jet resolution cut i/c. We 
see that for 0.04, Apo is constant within the errors. For smaller i/c there is a much 
stronger dependence. In this small yc region we may expect sizeable corrections involving 
Inycut- An attempt to resum leading and next-to-leading logarithms to all-orders for Apo 
for R2{D) will be discussed in sub-section 3D. In Figures 4-7 we used published OPAL 
data PI without correcting for hadronization effects. These are expected to be totally 
negligible for y^. ^ 0.04, but will be important in the small yc region. For thrust and EEC 
the ApQ^^ distributions are flat away from T ^ 1 and x ~ 180° where large corrections 
involving ln(l — T) and ln(cos^ |) respectively will be important. 

Returning to Figure 3 we see that there is significant scatter, in particular between 
the jet rates with EO and D recombination algorithms. To proceed further then we need 
either NNLO calculations or information on the Q-dependence of these observables. The 
only NNLO calculation completed so far is that for Rz [0. For i?3(£^0) and EEC there 
are published data from the JADE collaboration at the PETRA e'^e~ machine with ^/s 



in the range 22-44 GeV [^2|, using a comparable analysis. In the next sub-section we 
shall attempt to use this data together with the LEP data for these quantities to obtain 
absolute estimates of Apo from the Q-dependence. 

We conclude by tabulating in Table 2 the A^^ values obtained by fitting to LEP data 
at NLO for these same observables at the same values of the kinematical parameters as 
in Table 1 using p = Mz and p = Pec- The p = Pec values correspond to the A)-f- 
values for which Apo = 0, and so the scatter is directly related to points in the Apo plot 
of Figure 3. The p = Mz values show much greater scatter, since as discussed in section 2 
the scatter in r^^^{Mz) is superimposed on top of the physically- relevant scatter in Apo. 

To further emphasise this, suppose that we lived in an 'asymptotic world' where Q = 
Mz was sufficiently large that Apo ~ for all observables. Thus for all observables 

F{Rc,p{Mz)) = Mn - rp. (3.28) 



Supposing that A-^ = 200 MeV, then ( |3.28|) may be used to generate data 'RexpiMz 



in such an 'asymptotic world'. By construction the A^^;^ obtained by fitting at NLO with 
fJ' = fJ'EC will be 200 MeV for all observables, but that obtained with p = Mz will still 
exhibit significant, but completely predictable, scatter. These 'asymptotic world' results 
are also tabulated in Table 2. 



C. Apo from NNLO calculations and Q-dependence. 
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If we have available a NNLO perturbative calculation then the RS invariant p2 defined 
in equation (2.28) can be obtained. For i?^, p2 = —15.1 {Nf = 5) using the NNLO 
calculation of reference |T^. For all the other observables NNLO corrections have yet 



to be calculated and so p2 is unknown. If we insert the perturbative expansion of p{x) 
(equation (2.26)) into (3.8) we obtain 



R 



dx- 



(P2 + p3X + ...) 



P2R + 



1 + CX){1 + CX + P2X^ + . . .) 

• 2 



cp2)i?2 + 0(i?3)^Apr^° + 0(/?2 



(3.29) 



Using the 1993 LEP average data we have Rz = (19.97 ± 0.03)(1 + 6qcd) with 6qcd = 
0.040 ± 0.004 0. Thus with p2 = -15.1 we have Ap^^^^ ~ P25qcd ^ -0.60 ± 0.06. 
Adjusting the reference value of so that Apg^^ for Rz corresponds with this Ap^^^*^ 

results in A^^ = 288 ± 200 MeV. This is of course just the A;^ obtained at NNLO using 
the EC scheme with r = po and C2 = P2- 

To obtain estimates of Apo for the other observables one can try to use their Q- 
dependence. Suppose that we have measurements of R at two energies Q = Qi and 
Q = Q2 {Qi > Q2), then we can construct Apo^^(Qi) and Apo^^(Q2) of equation (3.27) 
from the data (with the same reference value of A^^). Then 



Ap^^^(QO - Apr (Q2 



bln^-F{R{Q^)) + F{R{Q2)) 



Q2 
R{Qi) 

R{Q2) 



I{x)dx = (RiQi) - RiQ2))IiR), (3.30) 



where 



I(x) 



p{x] 



+ 



1 



(3.31) 



is the integrand of Apo and R{Qi) < R < R{Q2 
measured from the Q-dependence of the data, 



x^(l + cx) 

Thus the integrand of Apo may be 



Apr(Q2 



iRiQ^) - RiQ2)) 



(3.32) 



Obviously by measuring dR/d\nQ a.t Q = Qq with sufficient accuracy one can in princi- 
ple determine I{R{Qo)) and so by suitably detailed measurement over the energy range 
[Q15Q2] one can determine I{x) on the interval [R{Qi), R{Q2)]. 

The situation for measurements at Qi and Q2 is shown in Figure 8 for I{x) versus 
X. The uncertainty in R is represented by the horizontal error bar between R{Qi) and 
R{Q2) and the measurement errors in R{Qi) and R{Q2) themselves contribute a vertical 
error bar. 
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We have from equation (3.29) that 

J(x) =p2 + (P3-2cp2)x + 0(a;2). (3.33) 

So /(O) = p2- Thus from a NNLO perturbative calculation we can obtain the integrand 
at the origin. Notice that, like A, the RS invariant p2 is connected with the asymptotic 
Q-dependence of R{Q). A next-NNLO calculation would provide ps and tell us the slope 
of I{x) at the origin, ps — 2cp2. 

If NNLO perturbation theory is adequate to determine Apo, then Apo ~ Ap^^-^*^ = 
P2-R(Q), corresponding to I{x) ^ p2 = constant, on the range [0,-R(Q)]. Thus, if NNLO 
perturbation theory is adequate, then we expect I{R) ~ P2 and so from the measurements 
of R{Qi) and /^(Qa) we can estimate Ap;^^*(gi) = I(R)R{Qi) from (3.32), 

A,est(n \ - (Apo"^(gi) - Apg"^(g2)) ^ , 

If NNLO calculations eventually become available for the observables one can then 
check explicitly whether Ap^^^'^ ^ Apg**. A marked discrepancy would indicate the 
importance of next-NNLO and higher perturbative effects and/or large 'non-perturbative' 
g-i/i? gffgc^s_ Given p2 one can then estimate pa from the slope of the straight line joining 
J(0) and I{R) and obtain an improved estimate of Apq'^*(Qi) from the area under the 
trapezium (see Figure 8). 

In this way perturbative calculations and experimental Q-dependence measurements 
serve as complementary pieces of information. The NNLO and higher perturbative cal- 
culations effectively provide details of the Q ^ oo running of R{Q) which could never 
be obtained experimentally. Together they can help to constrain the behaviour of the 
function ^{R) = —bp{R), and hence to refine the estimates of A^^. 

The running of the observables with energy has been used before as a test of QCD. 
For instance R3{E0,yc = 0.08) has been studied as a function of Q over the PETRA- 
LEP energy range [0. The apparent running has been compared with the NLO QCD 
expectations for various choices of renormalization scale p, and A. It has even been 
suggested that the LEP measurements may indicate slightly less running than is 
expected from QCD, and that light gluinos, which would modify the QCD /3-function, 
cannot be excluded. Such a claim seems ludicrously premature given our lack of knowledge 
of NNLO and higher QCD effects. 



If we take data for R^iEQ, = 0.08) at Q2 = 34 GeV from JADE |2[, and at Qi = 91 
GeV from OPAL [Q, then we obtain from (3.34), ApQ**(91) = — 5.5±3, which corresponds 
to I{R) = — 110 ± 60. So if NNLO perturbation theory is adequate for this observable we 
estimate p2 = —110 ± 60, for the as yet uncalculated NNLO RS invariant. For the EEC 
with X = 60° measured by JADE at Q2 = 34 GeV [||], and by OPAL at Qi = 91 GeV 



| 38| we similarly find ApQ**(91) = —1.5 ±1.5, which corresponds to I{R) = —21 ±21, and 



an estimate of p2 = —21 ± 21. Since ApQ^^ is very insensitive to yc and X; comparable 
results are obtained for other choices of these variables. 
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In figure 9 we plot Apo for Rz, R3{E0,yc = 0.08) and Reec{x = 60°) with Q = Mz- 
The diamonds correspond to absolute predictions. For Rz we take Ap^^-^"^ = —0.60 ± 
0.06, and RsiEO^y^ = 0.08) and for Reec{x = 60°) the Apg"*(91) values noted above 
obtained from Q-dependence. The crosses correspond to the Apg^^ for these observables 
with the reference chosen as 287 MeV so that Apg^^ for Rz agrees with Ap^^^*^. 
There is then agreement, within the considerable errors, between Ap^^^ and Apg** for 
i?3 and Reec- Demanding consistency between Apg^^ and Apg**, Ap^^^^^ then requires 
Ag) = 287 ± 100 MeV. 

Although only performed for a limited number of observables, we regard this con- 
sistency between Apo^^ and the absolute Apo obtained from NNLO calculations and 
Q-dependence measurements as very encouraging. Ideally one would like to see addi- 
tional LEP measurements for all observables off the Z peak at a lower value of Q to avoid 
the uncertainties due to combining measurements from different machines and detectors. 
One could then perform an exhaustive analysis of the kind described here and obtain a 
reliable determination of A^^. Scatter in the A^;^ values obtained by adjusting ApT^ to 
agree with Apg^* (obtained from Q-dependence using (3.34)) would then unambiguously 
indicate the importance of next-NNLO corrections or non-perturbative e~^/^ effects for 
some of the observables. 

It is interesting that the above Q-dependence estimates of p2, and the exact NNLO 
calculation of p2 for Rz, suggest p2 large (0(10)) and negative. This has implications for 
the infra-red behaviour of R{Q) since it is consistent with fixed-point behaviour, that is 



p{R*) = implying R ^ R* as Q ^ |31, 45 



D. Apo from next-to-leading logarithm resummations 



There has been much recent interest, both theoretical |l20| , pi[] and experimental in 
the possibility of resumming leading and next-to-leading logarithms (NLL) in kinematical 
variables to all-orders for LEP observables. By identifying the leading and next-to-leading 
logarithms in the perturbative coefficients = A^L'^^ + B^L'^^'^ + . . ., with L = In^ 
(where A is ?/c, 1 — T, cos^ |, . . .), and by demonstrating generalised exponentiation one 
can resum these terms to all-orders. To include the exact NLO perturbative coefficient, 
however, some ad hoc matching prescription is then required 0, and anyway the 
problem of renormalization scale dependence still remains. 

The advantage of the EC formalism here is that from the leading and next-to-leading 
logarithms in Tk one can obtain the leading and next-to-leading logarithms in the RS 
invariants pk = AkL?'^ + BkL"^^'^ + 0(L^^^^). One can then unambiguously construct 

p^L\a:) = a;2(l + ca; + p^^V + ... + p^^V...), 

pf^^ = AkL^" + BkL^"-'. (3.35) 
pNLL/^\ -g constructed from RS-invariants and so does not involve the renormalization 
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scale, or the NLO perturbative coefficient. There is, therefore, no analogue of the matching 
prescription ambiguity and one can obtain 



rRexpiQ) 

Ap^^^= / dx 
Jo 



^NLLi 



x) ~'~ x'^{l + cx) 



(3.36) 



This may be directly compared with ApQ^^. If the observed A-dependence of ApQ^^ in 
the small-A region disagrees with that predicted from Ap^^^ then one has unambiguous 
evidence that the NLL approximation is inadequate, that is that the neglected sub- leading 
logarithms and constants are important. 

We shall perform this comparison here only for the 2-jet rate i?2(-D), where the expo- 
nentiation of logarithms is rather straightforward, deferring discussions of other observ- 
ables. We have 

R2iD) = l + roR (3.37) 
where R has a perturbation series of the form (2.1). To NLL accuracy |20|, 



R2{D) = exp 



3CfL 



1 a 

3. 



6 



(3.38) 



where a = a{p = Q),L = \n^. The b dependent term in the exponent may be absorbed 

into a by a change of scale Q y\^^Q- Since the RS invariants pk are independent of the 
renormalization scale to construct them, we can equally use 



R2{D) 

One has from (3.37) that 



exp 



3C^L 



= exp(r, 



NLL^ 
^ 



R 



R2{D) - 1 



and from (2.27) with (3{a) = a^{l + ca) 



p(R) = a\R)(l + ca(R))^. 

da 



(3.39) 



(3.40) 



(3.41) 



From the form of pk of (2.28) one sees that C2, C3, . . . will not contribute at NLL level. We 



NLL I 



retain c in /3(a) so that p 



then using (3.41) one finds 



p^^^{R) = a'{R){l + ca{R))exp{r'^'^'^a{R)) 



X 



X {1 + cx + . . .). From (3.39) and (3.40) one has 



dR 

da 



exp^r. 



NLL^ 
^ 



NLL^ 



(3.42) 



(3.43) 
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where 

. Mlljl^. (3.44) 

If one then inserts p^^^{x) of (3.43) into (3.36) and change variables from x to 

ln(l + rf^ 

Performing the integration one finds 

NLL 

Ap^^ = F{u{Rexp)) - F{R,,p) - (3.46) 

In Figure 10(a) the diamond points show the absolute prediction for Ap^^^ versus ycut 
obtained using (3.46) together with the OPAL data of reference for R2{D) uncorrected 
for hadronization effects. The other points show the corresponding Apo^^ (see Figure 
5) with the reference value of adjusted so that Apo^^ agrees with Ap^^^ at the 
peak in Ap^^^. Notice that for yc ~ 0.01 the NLL prediction has a completely different 
low-j/c behaviour from the uncorrected data. In Figure 10(b) the same quantites using 
hadronization-corrected OPAL data for R2{D) |^ are displayed. Evidently, the low yc 
behaviour is now at least in qualitative agreement. We conclude that in order to test 
the adequacy of the NLL approximation we are having to compare with data in a region 
where hadronization corrections are sizeable. Unless one is extremely confident in one's 
ability to estimate these corrections it is hard to draw firm conclusions. 
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4 Summary and Conclusions 



Evidently if one has a NLO perturbative QCD calculation available there will be some 
uncertainty in the NLO perturbative prediction for a QCD observable due to the missing 
uncalculated higher-order terms in the perturbation series. This has the effect that when 
one compares the NLO calculation with experimental data one may anticipate that the 
value of Ajg^ extracted will not be universal but will exhibit some scatter for different 
observables due to the differing sizes of the uncalculated contributions. The problem 
is compounded by the fact that the Aj^ extracted depends also on the choice of RS at 
NLO, which may be labelled by ri the NLO coefficient (or equivalently the renormahsation 
scale fi), an unphysical parameter. We have shown in equation (3.25) that the relation 
between the extracted A{NLO,ri) and the actual which one is trying to measure 
is factorized into two contributions. One, AF{ri,R), is ri-dependent and known exactly 
from equation (3.26), and the other, Apo, is unknown but ri-independent. Apo depends 

only on the observable and NNLO and higher RS- invariants P2,P3, The predictable 

ri-dependence of A{NLO,ri) therefore has nothing to tell us about the importance of 
uncalculated corrections, the irreducible uncertainty resides in the unknown Apo- Some 
of the ri-dcpendent logarithms in the NNLO and higher corrections can be summed 
up into AF{ri,K) and the remainder can be absorbed into RS-invariant combinations 
contained in Apo. 

The unknown Apo can be isolated by choosing the particular RS where ri = 0, fi ~ 
^EC the effective charge scale, which sets AF = 0. From the scatter in A(A^LO, 0) for 
different observables one can then infer the scatter in Apo- in particular relative differences 
Apoyi — Apos for observables A,B can be absolutely measured. Figure 3 indicates that 
these relative differences cannot be neglected for a range of LEP observables. For at least 
some of the observables Apo must be sizeable, and so Apo must be estimated before we 
can determine Ajjg with any reliability. 

In contrast, in the standard LEP determinations of as{Mz) one tries to estimate 
the importance of uncalculated higher-order corrections by using two different ad hoc 
scale choices and interpreting the spread in extracted as{Mz) as indicating a 'theoretical 
error'. By artificially enlarging uncertainties in this way one obtains a spurious consistency 
between different observables, the real scale-independent uncertainty due to Apo being 
buried beneath the supposedly informative scale dependence of AF. The global as{Mz) 
determinations obtained with these sort of analyses should therefore be treated with some 
scepticism. 

Having determined that Apo is not negligible one must try to estimate it to make 
further progress. Within the effective charge formalism one can write A;^^ in terms of 
the measured observable R{Q), r^^^ and Apo constructed non-pcrturbatively from the 
measured running dR/dlnQ = —bp{R). Apo may also be obtained perturbatively by 
expanding the effective charge /3-function, so that given a NNLO calculation one can 
estimate Ap^^^^ — P2R. By combining NNLO calculations and Q-dependence mea- 



29 



surements one can refine ones knowledge of Apo and hence of A.jjg. The perturbative 
RS-invariants po{A), p2, ps, ■ ■ ■ are connected with Q ^ oo evolution of the observable 
R{Q), and hence provide information about the function p{R) in the vicinity of i? = 0. 

The present situation is that a NNLO calculation is available only for the hadronic 
width of the By estimating Ap^^^^ for this observable, and combining PETRA 
data for jet rates and energy-energy correlations, with LEP data on these observables, to 
obtain an estimate of Apo from Q-dependence, we found consistency for = 287 ± 100 
McV. 

We also showed that for observables where leading and next-to-leading logarithms in 
kinematical variables can be resummed one can unambiguously estimate Ap^^^ by resum- 
ming leading and next-to-leading logarithms in the RS invariants pk- In the conventional 
approach one needs to use an ad hoc matching procedure to include the exact NLO per- 
turbative coefficient, and the problem of scale dependence still remains. Performing this 
analysis for i?2(-D) we found qualitative agreement between the 7/c-dependence at low yc 
of Ap^^^ and Apg^^ using hadronization corrected data, but a completely opposite low 
He dependence using uncorrected data. 

We conclude that reliable measurements of Ajg^ at e+e" machines will require at least 
NNLO perturbative calculations and/or measurements of observables at more than one 
energy. Acquiring such information will entail considerable experimental and theoretical 
effort. The effective charge formahsm allows one to efficiently harness this hard won 
information to refine one's knowledge of Aj^g and the interplay between perturbative and 
non-perturbative effects. The insistancc on formulating everything in terms of physical 
quantities allows one to quantify uncertainties in a way which is impossible in approaches 
which choose the unphysical RS-dependence parameters according to some plausiblity 
argument. 

An important remaining problem is to extend the approach to processes with initial 
state hadrons where there is an additional factorization ambiguity connected with the 
separation of structure functions from the hard cross-sections. 
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Figure Captions 



Figure 1. The generic approximately inverted-parabolic behaviour of the NLO approxi- 
mant R^{t) versus a(r). The figure is plotted for po = 25 and Nf = 5. The dashed 
line represents measured data (see text). 

Figure 2(a). as{Mz) obtained by fitting NLO perturbative calculations with fi = Mz 
to OPAL data. 

Figure 2(b). As (a) but with enlarged 'theoretical errors' using /i = Mz and fits to data, 
as described in the text. Reproduced from reference [8]. 

Figure 3. Apg^^ of equation (3.27) for various LEP observables (see text for details). A 
reference value of A-V-ttt =110 MeV is chosen. 

MS 

Figure 4. The Ucut dependence of /S.p^Q^ for i?2(-E'0). OPAL data of reference [2] is used 
uncorrected for hadronization. A reference value of A-^ = 110 MeV is chosen. 

Mb 

Figure 5. As Figure 4 but for R2{D). 
Figure 6. As Figure 4 but for R^{EQi). 
Figure 7. As Figure 4 but for R^(D). 
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Figure 8. Learning about the integrand of Apo, of equation (3.31). NNLO pertur- 
bative calculations can provide /(O) = p2, and measurements at energies Qi and Q2 
can provide a point away from the origin. 

Figure 9. Absolute estimates of Apo for LEP observables at Q = 91 GeV are com- 
pared with the observed Apo^^ with = 288 MeV. For EEC(x = 60") and 
Rs{EO,yc = 0.08) the absolute estimates arc Apg''* of equation (3.34) obtained from 
Q-dependence using JADE data. For Rz Ap^^-^*^ from the NNLO calculation of 
reference [12] is used. 

Figure 10(a). For the observable i?2(-D) the NLL approximation Ap^^^ of equation 
(3.46) versus t/c is compared with ApQ^^ with A^^ adjusted. OPAL data of refer- 
vence [2] uncorrected for hadronization effects is used. 

Figure 10(b). As (a) but with hadronization-corrected data. 
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Obsorvril)lo 




,.:1/.S7 1 r \ 


T 


T— n Q 




R2{Eq) 


Vcut = 0.08 


10.22 




Vcut = 0.08 


10.09 


R2{D) 


Vcut = 0.08 


7.23 




Vcut = 0.08 


7.15 


EEC 


X = 60° 


11.30 


AEEC 


X = 60° 


5.13 


Rz 




1.411 



Table 1: NLO coefficients r^^{Mz) from [11] for various LEP observables at tfie partic- 
ular values of the kinematical variables noted. 



Observable 


LEP 

MS 


LEP 

MS 
1^ = I^EC 


Apo~0 
= 200 MeV 
^i = Mz 


Apo ~0 
Ag^ = 200 MeV 

= fJ^EC 


T 


883.0 


197.9 


894.4 


200 


R2{.Eq) 


335.7 


140.4 


498.5 


200 


RsiEo) 


326.8 


139.5 


488.3 


200 


R2{D) 


420.3 


252.1 


328.3 


200 


RsiD) 


425.6 


257.4 


325.2 


200 


EEC 


839.1 


268.2 


597.4 


200 


AEEC 


89.23 


63.79 


294.2 


200 


Rz 


254.4 


245.7 


206.9 


200 



Table 2: The first two columns show the central values of A^^ (in MeV) extracted 
by comparing NLO perturbative QCD predictions for the observables of Table 1, taking 
H = Mz and ^ = ^ec^ with OPAL data [2]. The third and fourth columns give the NLO 
A^^ values extracted with ji = Mz and jj = HeCi assuming Apo = for each observable, 
and an actual value of A^ = 200 MeV. 
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